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Abstract 

We present preliminary results of a search for B — > Xdj decays with a hadronic mass 1.0 GeV/c^ < 
M{X(i) < 1.8 GeV/c^. We consider seven final states with up to four charged pions and one 
neutral pion or r], which correspond to about 50% of the total fragmentation in this mass 
range. Based on a sample of 383 million BB events collected by the BaBar experiment at PEP- 
II, we measure a partial branching fraction J^l^^^^BiB ^d7)l(i.0GeV/c2<M(Xd)<i.8GeV/c2) = 
(3.1 ± 0.9^0^5 ± 0.5) • 10~^, where the uncertainties are statistical, systematic and model-dependent 
respectively. 
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1 INTRODUCTION 



We present an experimental study of the rare decays B X^j- Within the standard model 
of particle physics (SM), these flavor-changing-neutral-current (FCNC) b ^ d'j transitions are 
forbidden at tree level; the leading-order processes are one-loop electroweak penguin diagrams (see 
Figure [1]), where the top quark is the dominant virtual quark contribution. In the context of theories 



b 




d d 

Figure 1: Feynman diagram for ab^d'j transition 

beyond the SM, new virtual particles may appear in the loop, which could lead to measurable effects 
on experimental observables such as branching fractions and CP asymmetries [T]. 

Previous measurements of the exclusive decays B~^ P^l-, B^ p^'y, and B^ — > ujj by the 
Belle [2] and BABAR [3] experiments found branching fractions of 0{W~^), in good agreement with 
SM predictions [5]. In combination with the well-measured B — > K*'y branching fractions they 
also yielded measurements of the ratio of Cabbibo-Kobayashi-Maskawa (CKM) matrix elements 
iVfd/Vfsl which are in good agreement with results independently obtained from the ratio of -B^ 
and Bg mixing frequencies [1]. 

The experimental and theoretical uncertainties associated with the exclusive measurements are 
still large and there is a strong motivation to extend the experimental study to additional final 
states and different regions of the hadronic mass spectrum, with a measurement of the inclusive 
6 — > (i7 transition rate as the ultimate goal. In this note we report the first study of 6 — > (i7 
transitions using a sum of seven exclusive X^'j final states in the previously unaccessed hadronic 
mass range l.OGeV/c^ < M{Xd) < l.SGeV/c^. 

2 THE BABAR DETECTOR 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy e"*"e~ storage ring. Charged particle trajectories are measured using a five-layer silicon 
vertex tracker and a 40-layer drift chamber in a 1.5-T magnetic field. Photons and electrons 
are detected in a CsI(Tl) crystal electromagnetic calorimeter (EMC) with photon energy resolution 
aE/E = 0.023(£^/GeV)-i/'^e0.019. A ring-imaging Cherenkov detector (DIRC) is used for charged- 
particle identification. In order to identify muons, the magnetic fiux return is instrumented with 
resistive plate chambers and limited streamer tubes. A detailed description of the detector can be 
found elsewhere [6j. 
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3 EVENT RECONSTRUCTION AND SELECTION 



We exclusively reconstruct seven b ^ dj decay modes with up to four charged pions and one neutral 
pion or eta in the final state: 7r^7r~7, — > 7r"'"7r''7, B~^ — > 7r^7r~7r"'"7, B^ — > 7r^-7r~7r'^7, 

B^ — > 7r"'"7r~7r"'"7r~7, B^ 7r"^7r^7r+7r'^7, B~^ — > tt^7]j |8]. 

An important source of background is due to continuum events {e~^e~ qq, with q = u,d,s,c) 
that contain a high-energy photon from vr'^ or rj decays or from initial-state radiation. There is 
also background from combinatorial B decays. These include B X^tt^ and B X^r] processes 
that produce a high-energy photon in the tt^ or ry decay, as well as S — > Xgj decays where a is 
misidentified as a tt^. Background suppression cuts have been optimized for maximum statistical 
sensitivity S/y/lS + B), where S and B are the rates for signal and background events, using 
Monte Carlo (MC) simulated event samples, and assuming an inclusive 6 — > ^7 branching fraction 
of 1.0 X 10~^. Of particular concern are backgrounds from B meson decays that resemble the signal. 

The photon from the signal B decay is identified as a localized energy deposit in the calorimeter 
with energy 1.15 < E* < 3.5 GeV in the center-of-mass (CM) frame. The energy deposit must 
not be associated with any reconstructed charged track, be well isolated from other EMC deposits, 
and meet a number of additional requirements designed to eliminate background from hadronic 
showers and small-angle photon pairs. We veto any photons that can form a vr*^ (77) candidate by 
association with another detected photon of energy greater than 30 (250) MeV by requiring that 
the two-photon invariant mass be outside the range of 105 MeV/c2-155 MeV/c^ (500 MeV/c^- 
590 MeV/c2). 

Charged pion candidates are selected from well-reconstructed tracks with a minimum momen- 
tum in the laboratory frame of 300MeV/c. In order to reduce backgrounds from charged kaons 
produced in 6 — > 57 processes, a vr^ selection algorithm is applied, combining DIRC information 
with the energy loss measured in the tracking system. At a typical pion energy of 1 GeV the pion 
selection efficiency is over 85% and the Kaon mis-identification rate is less than 2%. 

vr^ [t]) candidates are formed from pairs of photons with energies greater than 20 MeV, with 
an invariant mass 117 < m^^ < 145 MeV/c^ (470 < m^^ < 620 MeV/c^). We require vr*^ and r] 
candidates to have a momentum greater than 300 MeV/ c. 

The selected tt^, vr", r/, and high-energy photon candidates are combined to form B meson 
candidates consistent with one of the seven signal decays. Here, the charged pions are combined 
to form a common vertex, where we require the vertex probability be greater than 2%. The mass 
of the hadronic system is required to be in the range 1.0 GeV/c^ < M{Xii) < 1.8 GeV/c^. We 
define AE = E*^ — -Ebeam' where E*^ is the CM energy of the i?-meson candidate and -B^eam 
is the CM beam energy. We also use the beam-energy-substituted mass ttt-es = \J -E'boam ~ Pb^ ' 
where p^ is the CM momentum of the B candidate. Signal events are expected to have a AE 
distribution centered at zero with a resolution of about 30 MeV, and an rriEs distribution centered 
at the mass of the B meson, m^, with a resolution of 3 MeV/c'^. We consider candidates in the 
ranges —0.3 GeV < AE < 0.2 GeV and ttt-es > 5.22 GeV/c^ to incorporate sidebands that allow 
the combinatorial background yields to be extracted by a fit to the data. 

Contributions from continuum background processes are reduced by considering only events for 
which the ratio R2 of second-to-zeroth order Fox- Wolfram moments [7] is less than 0.9. To further 
discriminate between the jet-like continuum background and the more spherically-symmetric signal 
events, we compute the angle 9t between the photon and the thrust axis of the rest of the event 
(ROE) in the CM frame and require \cos{Ot)\ < 0.8. The ROE is defined as all the charged tracks 
and neutral energy deposits in the calorimeter that are not used to reconstruct the B candidate. 
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The quantity cos{9j') and twelve other variables that distinguish between signal and continuum 
events are combined in a neural network (NN): R2, which is R2 in the frame recoiling against the 
photon momentum, is used to suppress events with initial-state radiation. We also compute the 
5-meson production angle 0|j with respect to the beam axis in the CM frame, and the Legendre 
moments Li = J2jPj ' \ cos 9j\^/J2jPj and Li = J2jPj ' I sin0||VX]jPj) where p* and 9* are the 
momentum and angle with respect to a given axis, respectively, for each particle j in the ROE. 
We use Li with respect to the ROE thrust axis, as well as L2, L^, L2, and L3 with respect to 
the high-energy photon direction as NN input. Differences in lepton and kaon production between 
background and B decays are exploited by including several flavor-tagging variables described in 
|10) . We select events for which the NN yields an output value greater than 0.83; this corresponds 
to signal and continuum background efficiencies of about 50% and 0.5% respectively. 

The expected average number of candidates per selected signal event is 1.75. In events with 
multiple candidates the one with the reconstructed TT^irj) mass closest to the nominal [9j is retained. 
For events containing no TT^{r]) the candidate with the highest vertex probability is retained. 

From simulated signal events with l.OGeV/c^ < Mx^ < 1.8GeV/c'^, applying all the selection 
criteria described above, we find an overall signal selection efficiency of 7.5%. 



4 MAXIMUM LIKELIHOOD FIT 



The signal content of the data is determined by means of a two-dimensional unbinned maximum 
likelihood fit to the AE and mES distributions. We consider the following contributions: signal, 
combinatorial backgrounds from continuum and BB background processes, B — > X^tt^ /rj decays, 
S — > decays, and self-crossfeed from mis-reconstructed B — > X^'y decays. The likelihood 
function is defined as 



exp 



hyp 



i=l 



n X! niViixj; Si 
j=i \i=i 



(1) 



Here, are the event yields for each of the signal and background hypotheses and N^at is the 
number of events observed in data; Pi is the probability density function (PDF) for hypothesis i, 
which depends on the fit observables xj = (A£^,mEs) and a set of parameters cJj. 

The functional form of each PDF is determined from MC simulated events. For the B — > X^'j 
self-crossfeed component we use a binned two-dimensional (Ai?, ttt-es) distribution as the input 
PDF, while for the other hypotheses we obtain the PDFs from unbinned one-dimensional fits to 
the AE and ttt-es distributions respectively. For the signal, the mES spectrum is described by a 
Crystal Ball function [12], and the AE distribution is parametrized as an asymmetric, variable- 
width Gaussian f{x) = exp —{x — p)"^ /{2a\ + aL^ji{x — fi)'^) , where is the peak position of 
the distribution, aL,R are the widths left and right of the peak, and ol^r are measures of the 
tails on the left and right sides of the peak. The contributions from B —>■ X^tt^ /t] and B Xg'y 
decays are both modeled by adding a Gaussian distribution to an ARGUS function [11] for tties 
and a Gaussian distribution to the second-order polynomial describing the background for AE. 
The Gaussian peaks in AE are displaced by ~ BOMeV for B — > Xdir^/rj by the missing photon, 
and by ~ 100 MeV for B — >■ Xg'j by the kaon mis-identification. A combined PDF is defined for the 
remaining B backgrounds and continuum processes using the sum of a Gaussian and an ARGUS 
function for mEs and a second-order polynomial for AE. The Gaussian component allows for a 
small component of B decays that peaks in ttt-es but not in AE. 
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In the likelihood fit the m-Es ARGUS slope parameter and the AE polynomial coefficients of 
the combined background PDF are free parameters, likewise the yields for this background and for 
the signal. All other fit parameters are fixed to the values obtained from the MC samples. 

5 VALIDATION OF THE ANALYSIS PROCEDURE 

In order to validate the analysis procedure, we embed signal MC events in backgrounds randomly 
generated from the background PDFs. The embedded events include both the correctly recon- 
structed and the self-crossfeed signal contributions. For different numbers of embedded signal 
events we determine the pulls on the yield and the statistical error from the fit. We find a bias 
in the fitted yield of +24% of the statistical error. We correct the final fit result for this bias and 
treat it as a systematic error. 

To check our analysis on data we select B K*'y decays, using the K* decay modes K*~^ — > 
K+7r° and K*^ K+Tr'. The selection is the same as for the two signal modes ir~^TT 7 

and B~^ tt^tt^"/, except that the vr^ identification requirements are replaced with a charged kaon 
selection, and the mass range is 0.6GeV/c2 < M{Xs) < l.OGeV/c^. In the fit procedure described 
in the previous section we include a i? — > Xg'y self-crossfeed component and a i? — > Xgir^ /r] 
background, but not a B ^ X^j misidentification component because this is expected to be 
negligible. B K*j has a large signal yield, so we use this fit to determine the signal shape in 
data by allowing the means and widths of the signal PDFs to vary. We find 1680 it 51 i? — > K*'y 
events in the BABAR data, where this error is purely statistical. This is in excellent agreement 
with the expectation of 1616 it 28 events based on the previously measured branching fraction of 
13{B K*-f) = (4.02 ± 0.33) x IQ-^ [T3|. 

As a further check on data we perform the signal analysis for the decays B^ 7r+7r^7, B'^ — > 
TT+TfO'^^ and B^ ■K~^TT~7r^^ in the hadronic mass range 0.6GeV/c^ < M{X(i) < l.OGeV/c^, which 
contains the p^,p^ and u> resonances. In this fit we include all the contributions in a similar way 
to the B Xd'y analysis in the higher mass region, and fix the signal PDF shape to the values 
obtained from the B K*^ fit. We obtain a combined B — > {p,uj)^ signal of 73 it 25(stat.) 
events corresponding to a statistical significance of 2.9o". This is consistent with the 66 it 26 events 
expected from the average of previous branching fraction measurements by BABAR and Belle of 
B{B — > {p,u!)'y) = (1.28 lb 0.21) x 10~^ ^14j. Note that we do not make use of the resonance mass or 
the decay helicity angle to discriminate against backgrounds, so we do not expect this measurement 
to be as precise as the exclusive B {p,uj)^ analyses. Figure [2] shows a comparison of the PDF 
component shapes (curves) to the data (points). 

6 RESULTS 

Finally we perform the full fit for the sum of all seven decay modes in the hadronic mass range 
l.OGeV/c^ < M{Xd) < 1.8 GeV/c^, again using the signal PDFs from B K*j. 

Figure [3] shows a comparison of the PDF component shapes (curves) to the data (points). For 
each plot, a cut is applied to the variable not plotted around the signal peak. 

Figure 0] shows a comparison of the PDF shapes (solid curves) to the data using the event- 
weighting technique described in Ref. [13] to subtract the fitted background. For each plot, we 
perform a fit excluding the variable being plotted and use the fitted yields and covariance matrix 
to determine the relative probability that an event is signal or background. The distribution is 
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Figure 2: Projections of tues with -0.1 GeV < AE < 0.05 GeV (left) and AE with 5.275 GeY/cP < 
"T-ES < 5.286 GeV/ (right) in the B — > (p,a;)7 fit to data (points with statistical errors). The 
curves represent the PDFs used in the fit with the combined PDF in sohd blue. The dashed 
lines represent the individual PDFs with with signal component in green, signal cross-feed in pink, 
h background in blue, B Xdir^/r] background in brown and background in orange. 



normalized to the yield for the given component and can be compared directly to the assumed 
PDF shape. We find good agreement between the data and the PDFs. 

We find a signal yield of 178 =b 53 events and a background of 27670 ± 173 events in the full 
fit. Taking only statistical uncertainties into account, this corresponds to a signal significance of 
3.4(7. We calculate the partial branching fraction summed over the seven modes in the mass region 
1.0 GeV/ < M{Xd) < 1.8GeV/c2 using the following formula: 

7 ^ 

X! ~^ ^dl)\{1.0GeV/c'2<M(Xa)<1.8GeV/c'2) = 7, ^^^S (2) 

where ngig is the number of fitted signal events, e is the signal selection efficiency, n^-^ is the 
number of BB pairs in the dataset, and the factor 2 is to account for the possibility of either B 
decaying into a signal mode. We obtain the following result: 

7 

B{B ^ Xd7)l(i.0GeV/c2<M(x,)<i.8GeV/c2) = [3-1 ± 0.9(stat.)] • 10-6 ; (3) 
The relevant systematic uncertainties are discussed below. 

7 SYSTEMATIC UNCERTAINTIES 

Table [T] gives an overview of the contributions to the systematic uncertainties. These are associated 
with the signal reconstruction efficiency, the modeling of the signal and BB background PDFs in 
the likelihood fit, and the modeling of the B X^'y spectrum and the fragmentation of the X^ 
system. 

The systematic errors affecting the signal efficiency includes uncertainties on tracking, particle 
identification, 7 and vr" reconstruction, the tt^ j-q veto, and the neural network selection. These 
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Figure 3: Projections of mES with -0.1 GeV < AE < 0.05 GeV (left) and AE with 5.275 GeY/cP < 
iTT-ES < 5.286 GeV/ (right) in the 6 — > (i7 fit to data. The curves represent the PDFs used in the 
fit with the combined PDF in sohd blue. The dashed lines represent the individual PDFs with 
with signal component in green, signal cross-feed in pink, 6 ^ 57 background in blue, B X^tt^ /rj 
background in orange and background in brown. 




Figure 4: Signal distributions of rriES (left) and AE (right) for B X^'y in the range 1.0 GeV/ < 
M{Xd) < 1.8 GeV/ with the background subtraction as described in the text. The curves represent 
the PDFs used in the fit, normalized to the fitted yield. 



uncertainties are evaluated using independent data and MC-simulated event samples. They give a 
combined systematic error of 6.5% when added in quadrature. 

To estimate the uncertainties related to the signal and background yields we vary the parameters 
of the PDFs that are fixed in the fit. For the signal PDF the means and widths are varied within 
the range allowed by the fit to the B — > K*j data. Other fixed PDF parameters are fluctuated 
within their statistical errors. The combination gives errors of ±13% on the signal yield. 

All relative and absolute normalizations of background components that are fixed in the fit are 
also varied. The absolute normalizations of the B — > X^^y cross feed, B —f Xs^, and B X^tt^ /r] 
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Table 1: Fractional systematic errors of the measured branching fraction. 



Source 


Relative uncertainty (%) 


Tracking efficiency 


1.7 


Charged-particle identification 


2.0 


Photon selection 


2.5 


vr*^ and r] reconstruction 


1.7 


vr*^ and r] veto 


1.0 


NN efficiency 


5.0 


PDF shapes 


13.0 


B background normalization 


+10.5 
-3.0 


Fit bias 


7.1 


B counting 


1.1 


Combined 


+19.3 
-16.5 


Signal model 


3.4 


Xd fragmentation model 


13.6 


Combined 


14.7 



components are changed by ±50%, ±20%, and ±100% respectively; and the relative contribution 
of the remaining B decays to the combined background component is varied by ±20%. The full size 
of the bias on the signal yield observed in MC experiments (see Section [5]), is taken as a systematic 
uncertainty of the fit procedure. 

There is a small 1.1% uncertainty on the overall normalization associated with the the total 
number of BB pairs in the underlying data sample. 

The impact of the assumed spectrum on the efficiency is studied by using two different sets 
of values for the kinetic parameters of the Kagan-Neubert model [15] (m6,/x^) = (4.65, —0.50) and 
(4.80,-0.35). These are consistent with fits of B ^ Xgj data p!6]. The default fragmentation 
model for the X^ is JETSET [I7j , known not to give a good description of the distribution of final 
states in B ^ Xg^ [H]. To estimate the uncertainty associated with the fragmentation model, 
we re-weight the relative contributions of the seven signal modes according to results obtained in 
|18] for the corresponding B — > Xg^ decays, and take the resulting change in the signal yield as 
a systematic error. Combining these two uncertainties, we assign a systematic error of 14.7/the 
modeling of the signal. 

8 SUMMARY 

In this paper we have demonstrated the feasibility of measuring B — > Xd^f decays with a semi- 
inclusive technique using a sum of seven exclusive final state in the mass range l.OGeV/c'^ < 
M{Xd) < 1.8GeV/c^ and presented the first evidence for b ^ dj transitions in this mass range. 
We find 178 ± 53 signal events. 

We measure a partial branching fraction ^5cd=i ^(^ ^ -'^d7)l(i.0GcV/c2<Af(Xd)<i.8GeV/c2) = 
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(3.1 zbO-O^Qg ±0.5) -10 ^ , where the uncertainties are statistical, systematic and model respectively. 
Including statistical and systematic uncertainties only, this corresponds to a signal significance of 
S.lcr. 

9 OUTLOOK 

In the future, the inclusive b —f dj transition rate can be measured. This will require further study 
to: 

• Extend the current measurement to the mass region 0.6 GeV/ < M(Xrf) < 1.0 GeV/c^, which 
has been used as a control sample in this analysis with only three decay modes reconstructed. 
We can measure the seven decay modes over both ranges. 

• Correct for the part of the fragmentation that is not reconstructed by the seven decay 
modes. Based on MC simulation and the JETSET model we expect the measured seven modes 
to account for about 50% of the the hadronic mass region 1.0 GeV/ < M{X(i) < 1.8 GeV/ c^. 

• Extrapolate to the full mass region. Assuming the shape of the hadronic mass spectrum is the 
same for b ^ dj as for 6 — > 57 decays [18], the mass region 0.6 GeV/ — 1.8 GeV/ comprises 
roughly 60% of the spectrum. 

Further, the extraction of the ratio |Vj(i/Vjs| can be performed by comparing the branching 
fractions for b ^ d'j and 6 ^ 57 over the experimentally accessible final states and hadronic mass 
range. By measuring this ratio over a larger hadronic mass range, and with a larger set of final 
states it should be possible to reduce the theoretical uncertainties compared to the ratio of B ^ 
and B — > K*j decays. 
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